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Abstract 

By coupling axions strongly to a hidden sector, the energy density 
in coherent axions may be converted to radiative degrees of freedom, 
alleviating the "axion energy crisis" . The strong coupling is achieved 
by mixing the axion and some other Goldstone boson through their 
kinetic energy terms, in a manner reminiscent of paraphoton models. 
Even with the strong coupling it proves difficult to relax the axion 
energy density through particle absorption, due to the derivative na- 
ture of Goldstone boson couplings and the effect of back reactions 
on the evolution of the axion number density. However, the distri- 
bution of other particle species in the hidden sector will be driven 
from equilibrium by the axion field oscillations. Restoration of ther- 
mal equilibrium results in energy being transferred from the axions to 
massless particles, where it can redshift harmlessly without causing 
any cosmological problems. 



1 Introduction 



Goldstone and pseudo-goldstone bosons play a large role in the speculations 
of particle physicists and cosmologists. Some oft-discussed examples are 
axions, 1 majorons, 2 familons, 3 and the non-abelian goldstone fields of texture 
models. 4 The salient feature of the dynamics of goldstone bosons is that they 
couple derivatively. If the associated global continuous symmetry is broken at 
a scale /, then the coupling always involves powers of f~ 1 d fl b, where b is the 
goldstone boson. The consequence of this is that if / is very large compared 
to the relevant masses and momenta, the goldstone boson decouples. This 
fact is exploited in invisible axion models, 5 where for f a > 10 10 GeV the 
axion's couplings are small enough for it to have evaded detection and avoided 
conflict with astrophysical observations. 6 This fact also lies behind the well- 
known "axion energy problem", 7 which is that if f a £ 10 12 GeV the energy 
in coherent oscillations of the axion field after the QCD phase transition in 
the early universe is unable to dissipate due to the axion's weak coupling and 
eventually overdoses the universe. 

In this letter we show that different goldstone fields can mix with each 
other through kinetic terms in the Lagrangian in a way reminiscent of photon- 
paraphoton mixing, 8 and that a goldstone boson can thus acquire effective 
couplings at low energy to certain fields which are much larger than the 
typical p/f. After demonstrating this point we show how it makes possi- 
ble scenarios in which the energy in coherent oscillations of the cosmological 
axion field can be dissipated and the upper limit on f a removed. 
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2 Mixing of goldstone bosons 

Consider a theory with a global symmetry U(l) a x U(l)b- Let Q and $ be 
two scalar fields with charges (1,0) and (0,1) under this symmetry. Suppose 
that = f a and |(<E>)| = fa, where fa is assumed to be much smaller than 
fa, and call the resulting goldstone modes of f2 and <3>, 'a' and '6' respectively. 
Then it is clear that couplings of a will be suppressed by pj fa and those of 
b by pj fa. 

It can happen that by integrating out some fields of mass M an effective 
higher dimension operator will result of the U{l) a x [/(^-invariant form 

A™ = ^n%n&d»<!> + h.c. (i) 

If = (f a +Q)e ta lf a and $ = (fa+^)e lb ^ b are substituted into this expression, 
the result contains the term 

C a , b ^e(d^ad%) (2) 

with 

e - ■ ( 3 ) 

Through this mixing of a and b, the 'axion' a can couple to particles which 
have U(l)b charges with effective strength ep/fa, which can be much greater 
than p/ f a (see Fig. 1). 

For example, a term such as in Eq. (1) could arise from integrating out a 
field 7], with charges (—1/2, —1/2), mass M, and coupling [\(il&r] 2 ) + h.c.}. 
For simplicity let r\ have positive mass-squared and no VEV. Then the loop 
shown in Fig. 2 will produce the following effective term (among others) for 



momenta less than M: 

A 2 



192tt 2 M 2 
and therefore, from Eq. (3), 

A 2 fafb 



(Stf + h.c.) (4) 



(5) 



96tt 2 M 2 

Now, given that the term \(VL<^rf) exists, it would require an artificial can- 
cellation for M 2 to be much less than Xf a fb, and hence for e to be much 
larger than A/(967r 2 ) ~ 1CT 3 . There is also an absolute limit that e < 1, 
since the kinetic terms of a and b have the form 



A« i)(^) 



and a wrong-sign metric would otherwise result. There is no reason, however, 
why e cannot be of order, though less than, unity. 

If the higher dimension term in Eq. (1) were a "hard" term all the way up 
to momenta of order /„, then the theory would become strongly coupled and 
physics uncalculable unless M 2 were > / 2 and hence e were less than about 
fb/ fa- The effective coupling of a to particles with U(l)f, charges through 
mixing then would be of order ep/fb < p/ fa- 111 other words, the coupling 
of a would not be enhanced by this mixing. If, on the other hand, the term 
in Eq. (1) arises through a loop such as that shown in Fig. 2 where the 
intermediate particles have mass M 2 ~ fafb-, then it softens for momenta 
above M and no problem arises in higher order. 
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3 Alleviating the axion energy crisis 



We wish to make use of the mixing of goldstone bosons to let axions couple 
strongly to some other particles and thus, somehow, allow the primordial 
energy density in axions to be dissipated. 

Let the field a be the axion that solves the strong CP problem. Then 
U(l) a has a QCD anomaly and a gets a mass of order A 2 QCD / f a from instan- 
ton effects. Several astrophysical arguments 6 suggest that f a > 10 10 GeV, 
otherwise the evolution of stars or their remnants would be drastically differ- 
ent than observed. Thus the coupling of the axion to ordinary matter, which 
is proportional to 1/ f a , is highly suppressed. 

The axion energy crisis arises as a result of the small coupling and the 
small mass of the axion. At temperatures above the QCD phase transition the 
axion mass is small enough that a ~ 0. As the universe cools, m a increases 
until the axion field becomes dynamic, i.e. m a > H, where H is the expansion 
rate. Barring a coincidence of alignment for the Peccei-Quinn angle, the 
initial value of the axion field is a ~ f a . The axion number density is then 
n a pa f%H. After that point the number of axions in a comoving volume is 
constant, while the axion mass increases from H to its zero temperature value 
m a0 . If the energy density of axions after the QCD phase transition exceeds a 
critical value, p a ^ cr ~ 1CT 7 Aq C < d , then the the universe would become matter 
dominated too early, i.e. at temperatures > 10 eV. The energy density in 
axions at T ~ A QC d is roughly p a (A QC D) ~ A. A QCD f a /M Ph which results in 
the constraint f a < 10 12 GeV. In non-inflationary models there are bound 
to be axionic strings - whose decay may 9 or may not 10 increase the axion 
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abundance 100-fold, lowering the cosmological bound to f a £ 10 10 GeV. 

There have been several approaches to solving the axion energy crisis. 
For example, entropy production 11 may dilute the energy density in axions 
relative to ordinary matter, or inflation may fix the initial value of the Peccei- 
Quinn angle to a suitably small value 12 . However, attempts to dissipate the 
axion energy typically fail. Direct axion absorption presents several chal- 
lenges. First, the goldstone boson nature of axions demands a factor of 1/ f a 
at each vertex. Second, since the axions to be absorbed are non-relativistic 
the derivative couplings demand an additional factor of m a at each vertex. 
Third, any process for single axion absorption must admit axion emission in 
the time reversed channel. This leads to a reduction in the dissipation rate 
by a factor of m a /T from a naive approach that would ignore the back reac- 
tions. Flynn and Randall 13 considered the non- derivative coupling of axions 
to mesons 14 , but the first and third arguments above were sufficient to leave 
axions undamped by several orders of magnitude. 

In the present case, the mixing of axions with some other goldstone boson 
allows the possibility of coupling axions strongly to some other non-standard 
matter. However, without some explicit symmetry breaking of the i7(l)&, 
the axions must couple derivatively. The combination of powers of m a from 
the couplings and from consideration of back reactions has caused all our 
attempts at axion absorption to fail. 

Nonetheless, there is a mechanism whereby axion oscillations can be 
damped, allowing for f a of order the grand unified scale M GUT , or even 
the Planck scale M Pt . Let the axion be mixed, as described above, with 
strength e with a goldstone boson b arising from the breaking of C/(l)&. The 
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axion coupling to '^-matter', as we will call it, will be vastly stronger than its 
coupling to ordinary matter. Our basic idea is then that the oscillating axion 
field causes the energies of 6-particles to change in energy. If a thermal bath 
of such 6-particles were present in the early universe when the axion field 
started to oscillate, then those oscillations would drive the 6-distribution out 
of thermal equilibrium. By adjusting the interaction rates of the 6-particles 
to a suitable value the axion energy can be dissipated into thermal energy of 
b particles. The final step is to note that the b goldstone boson is massless, 
so that all b energy can in principle redshift away as radiation. The axion 
energy is thus put into a harmless form which never dominates the universe. 

4 An explicit model 

To illustrate how this may happen let us calculate in a simple model in which 
the '6-matter' consists of the field <3> whose vacuum expectation value breaks 
U(l)b and whose phase is 6, and a scalar field x which has no VEV. Let 
the U(l)b charges of $ and x both be +1 and let them have a coupling 
(\g<S>* 2 X 2 + h.c). 

The terms of the Lagrangian that will be relevant for our analysis are 
C = \{d,af-\mla 2 + l -{d,bf + e{d,ad^b) + 

\\d»X? - \m 2 \x\ 2 ' \f {e~ 2tb/fb X 2 + h.c) . (6) 

The last term arises from expanding $ as (fb + ^)e tb ^ b in the $* 2 x 2 coupling 
and defining /i 2 = gf^. At this point we are free to choose the parameters 
of C. For purposes of illustration it is convenient to choose /& ~ m ~ Aqcd 
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and to choose \i somewhat less than m. These values are not absolutely 
necessary, but they clearly show the damping mechanism to work. 

The first task is to resolve the mixing of a and b. To bring the kinetic terms 
of a and b to canonical form we perform a non-orthogonal transformation 



a' = ay/l — e 2 

b' = b + ea . (7) 



This gives 



£ = \W - iml(l - e 2 )-V s + \(d,Vf 

+ \Kx? - \™'\x? - \»> (e 2 " A| -* + 7^V + k-c) • (8) 

At this point it is convenient to drop the primes and work with the canon- 
ically normalized fields. For further convenience we represent the coherent 
oscillation of the axion field by 



a(t) = a Vl — e 2 sm(m a t) . (9) 

Then the effective Lagrangian for x m t ne presence of axion oscillations 
becomes 

Ceff = \\d,X? - \m 2 \ X ? - -/ {e mt) X 2 + h.c.) (10) 

with 0(t) = (ea /f b )sm(m a t). 

It is important to note that even though a(t) is slowly varying (with 
frequency m a ), 6{t) is rapidly varying: q = 6{t) = e^ h2 -cos(m a t). Since 
aom a = p l J 2 on ~ A 2 QCD (f a /M p i) 1/2 (initially), q ~ eA QC D- For time scales 
short compared to m~ l we may treat q as being constant and write C as 

C ef f = \\d,x? - \m 2 \x? ~ -/ [e 2tq Y + h.c] . (11) 



The problem now is to find the evolution of the x field, a non-trivial task 
because of the explicit time dependence. This is done in the Appendix, but 
we give a simplified version here. First, if fi 2 = then for each momentum 

— * 

k there would be two degenerate eigenstates. One could choose to work 
in a basis with linearly polarized or with circularly polarized states. With 
degenerate states it is easy to show that the Hamiltonian is constant for 
any linear combination. Next, consider the case /i 2 7^ 0, but q — 0. Then 
the two eigenstates have masses (m 2 ± /i 2 ) 1//2 , and correspond to the linearly 
polarized states along the real and imaginary components of x\ *.e. they 
are A(x^) and B(x^), where x( x>x ) = A(x^) + iB(x^). Now, there is no 
explicit time dependence so the Hamiltonian is still constant. This can be 
seen explicitly by writing the Hamiltonian in terms of the linearly polarized 
states and observing that there is no mixing of different frequency states in 
the Hamiltonian. Finally, for non-zero q, the Hamiltonian is time dependent 
and the energy of the x field varies with time. For fi 2 small the evolution 
is well approximated by the fi 2 = eigenstates, but there is a qualitative 
difference between the linearly polarized and circularly polarized solutions. 

For the linearly polarized case, the potential in the A (or B) direction 
has a curvature that varies sinusoidally in time: m 2 ^ = m 2 ± /i 2 cos(2gt). 
Since we will assume q <C m, the variation in Tn 2 e ^ is adiabatic. In that 
case one expects each quantum to have energy \Jm 2 + k 2 + /i 2 cos(2gt) ~ 
Uq + |^cos(2gt), where c^o = \J fn 2 + k 2 . The number of quanta of oscillation 
is an adiabatic invariant. 

On the other hand, for the circularly polarized case, the state samples 
the potential at all phases of x- Since m ^> q, this averaging takes place in 
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a time small compared to the evolution of the potential and there is little 
time dependence in the energy. In the Appendix we show that, in fact, the 
true eigenmodes are approximately circularly polarized. There is no time 
dependence for pure states; however, in a thermal bath we do not expect pure 
states unless there is some conserved quantum number. The corresponding 
symmetry in our case is the b symmetry, which is broken, so that b charge 
is not conserved. We therefore expect all polarizations to play a role, and 
having identified the linearly polarized states as having an energy dependence 
driven by the axion field we look for dissipation of the axion energy through 
those modes. 

Time variation of the energy levels of the x fi e ld can lead to dissipation. 

— * 

Suppose the rate at which a linearly polarized \ particle of momentum k 

— * 

scatters is given by 'j(k). Then the occupation number f%, of that state will 
obey the differential equation 



h = 7(*) 



(12) 



where E(k, t) = cj + |^cos(2gt) and / is the equilibrium occupation number 



(Bose-Einstein). Treating ^ = as small, this has solution 

/ E (^/„(^)+/^)^^cos( M <H» 

where f' Q = df /dx and tan/5 = —(2q)/j. The occupation numbers thus 
'lag' by a phase (5 behind the equilibrium value as that oscillates in time. As 
7 — > oo, this lag goes to zero and f%(t) — > fo ^ E - ^ ; the system remains 
in equilibrium. As 7 — ■> 0, the interactions turn off and there is no time 
variation of f%(t). Maximum dissipation occurs for 7 = 2q. The average 
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power dissipated over a period (2q) *by the mode with momentum k is given 
by 



P = -(2q) [ Q dtf % (t)E(k,t) 
f , /u \ 7 (2g) 2 p 4 



T) + {2q) 2 AujIT 

dt sin(2gt + /3)cos(2gt) (14) 



or 



/ojA 7 (2g) 2 /x 4 



The total power dissipated per unit volume is then given by 

7(2g) 2 



_1 /■ d 3 k fuo\ 
F ~ 3 J (2tt)3 /o [t) 



I* 



(16) 



_7 2 + (2g) 2 _ 

The factor 1/3 comes from averaging over all possible polarizations of the x 
particles. 

Equating the dissipation power with the change in axion energy density 
P = p a , and using (q 2 ) = ((ep)/ 2 / f b cos(m a t)) 2 ) = ^ 2 jrpa, the energy in 
coherent axion oscillations evolves as 

where C(m,T) = J dxj^f^c 2 + x 2 ), with c = m/T and x = k/T. Let 
us now make the approximation that whatever dissipation occurs takes place 
in less than an expansion period - so that the functions 7 and C(m,T) may 
be approximated as constant. Then, taking the integral / dt ~ O.IM^/T 2 , 
the final axion energy density is 

Paje^ = Paji e~F r e (18) 
10 



This solution embodies a number of special cases, but we are mostly in- 
terested in the cases where we can treat (2e 2 p a )/(/ 2 7 2 ) as small. Then, 
the axion density decreases exponentially during one expansion time. Given 
that the maximum initial value of p a is approximately Aq CD , as long as 
(10 _3 e 2 /i 4 M p i)/(/ 2 7T 2 ) > 20, the axion density will be reduced to accept- 
able levels. 

We take the temperature to be about Aqcd- Dissipation is fastest if 
fb is small. On the other hand, if /& < T, the U(l)b symmetry will be 
unbroken and the axion will have no goldstone boson, b, to mix with, and this 
dissipation mechanism will not operate. So the greatest dissipation should 
exist for fb ~ T ~ Aqcd, and 

c2 (it-)' (it-) ^ 10-16 ' (19) 

On the other hand, the dissipation cuts off for q <C /i 2 /m. So the energy in 
the axion field cannot be reduced below 

Pf,min « ( 20 ) 

by this mechanism. Demanding that this be less than ~ 10~ s Aq CD in order 
to comfortably solve the axion energy problem gives (/i/ 'Aqcd) 4 < e 2 lCT 8 , 
and from Eq. (19), e 4 (7 / Aqcd)^ 1 > 1CT 8 . These two constraints can easily 
be satisfied, for example by e ~ 10" 3 , 7 < 10~ 4 Aq C d, and fi 2 10~ 7 Aq CD . 

5 Discussion 

We have established that it is possible to dissipate the primordial energy 
density in axions through their interactions with other particles. A key fea- 
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ture of the process we envision is to mix the axion with some other goldstone 
boson b, whose decay constant /& is much smaller than f a . The mixing may 
be large, and this allow axions to couple strongly to L b matter'. By adjusting 
the matter content of the 6-sector we can arrange for the axion oscillations 
to drive oscillations in the energy of b matter. This pushes the b matter out 
of thermal equilibrium, and ultimately allows the axion energy density to 
dissipate in the 6-sector. 

Although, we have constructed a model whereby the axions dissipate, the 
addition of new particles may have other troublesome consequences, chief of 
which is that b matter may itself come to dominate the universe. The easiest 
way to avoid this is to adjust the masses and couplings of the 6-sector so the 
real part of the $ field is the lightest 6-particle other than the b goldstone 
boson itself. Then as the universe cools, first all b matter will end up as p&'s 
which will then decay with a rate of order into b goldstone bosons. 

We argued earlier that it is difficult to dissipate axion energy through the 
process of single axion absorption, yet, in the present scenario axion energy 
is indeed absorbed. It is natural to ask if there is a "Feynman diagram" 
explanation of our process. The answer is yes. Imagine a process whereby a 
single axion is absorbed from the coherent state, as in Fig. (3a) and a second 
process whereby two axions are absorbed, as in Fig. (3b). The ratio of the 
rates for these two processes is 

r = r 2 /r! ~ n a — / , (21) 
Jb 

where n a is the occupation number for the state. As stated earlier the effec- 
tive single axion absorption rate is reduced by a factor of m a /T when back 
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reactions are included. However, this is also true of the two axion absorption 
rate, and so the ratio in Eq. (21) is correct even in the presence of backre- 
actions. Now, for the coherent state, n a is approximated by n a ~ f 2 /H 2 , so 
the ratio becomes R ~ e 2 M 2 J Aq CD , which is quite large. This leads to the 
situtation that two axion absorption formally exceeds one axion absorption. 
The three axion rate would be even larger... In this situation what one must 
do is sum over diagrams with any number of axions attached to the particles 
participating in the scattering process. This is equivalent to solving for the 
propagator of these particles in the presence of the coherent oscillating axion 
field, as we have done in this paper. 

Along these lines, even in a normal axion model with no mixing each 
additional axion brings a factor of M v \ / f a to the amplitude, so here also one 
should not calculate individual absorption rates. Rather, one should calculate 
the propagators for the other particles, including the time dependent axion 
field, to arrive at the axion dissipation rate. We are presently considering the 
question of what dissipation may result from the non-derivative couplings to 
mesons, and plan to present results in a separate paper. Although derivative 
couplings may seem important too, one must remember that such couplings 
are less effective at damping the non-relativistic coherent axion oscillations 
due to the factor of m a in the coupling. 

Besides damping the coherent oscillations in the early universe, the kind 
of mixing of axions with other goldstone bosons we have been discussing 
could lead to other interesting effects. In particular, it is possible now to 
contemplate that even invisible axions with f a ~ M GUT or M Pi can have 
sizable couplings to some kinds of matter, perhaps even to particles that 
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carry standard model gauge charges. This would mean that the "invisible 
axion" may not be quite as invisible as was thought. 

Also worth further investigation is whether similar mixing in familiar 
majoron, familon or texture models could lead to interesting phenomena. 

Appendix 

The easiest way to find the eigenstates of the Lagrangian in Eq. (11) is to 
make the field redefinition if> = e iqt x, after which the Lagrangian becomes 

£' = i(|^| 2 -(m 2 -g 2 )|^| 2 + ^*-^-^ 2 (^ 2 + ^ 2 )) • (Al) 
The equations of motion for £ are 

i> - vV + ( m2 - o 2 )^ - 2i# + ifip* = o 

4>* - V V + (m 2 - q 2 W + 2 *#* + ffip = . (A2) 
The solutions to these equations are 



or 

= a,- 



where the index % denotes one of four solutions for a given wave vector k. 
Due to the mixing of ip and if}*, to simultaneously solve both equations of 
motion requires that ctj be either pure real or pure imaginary, but otherwise 
the normalization amplitude is arbitrary. The positive eigenfrequencies are 
given by 



U}± 



ul + q 2 ± 4u£q z + n 



.22 , ,.4 



1/ 



/21 V2 



(A4) 
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Negative frequency solutions for momentum k are equivalent to positive fre- 

— * 

quency solutions with momentum —k, and so are not distinct. For the case 
when ctj is real, the mixings are given by 

,2 



5 + = 



-/r 



For the case where ctj is imaginary, the sign of S is changed. 



(A5) 



Counting ou±, and the real and imaginary a iy there are eight eigenfunc- 
tions combining the wavenumbers k and —k, which is the right number to 
determine the field and its first time derivative at each point; so this is a 
complete set of states for the ip field. 

For the x held, the states are given by Xi = e~ tqt ipi- We can choose any 
four of the eight eigenstates to characterize the momentum k (using the other 
four for —k). It seems natural to chose a set where u — q — > oo as fi 2 — > 0, 
since in that limit we must get back the free field theory for x- With this in 
mind, the general solution that we will associate with momentum k is 



X=J2 ( a s+Xs+ + ia s -Xs~) 

s=± 

where the a,,/ are real numbers and where 



Xss > = N; s } /2 \(1 - 5 sg/ ) e *(«.-9)*-*« + 5 ss/e - i ^ +q)t+ik - x 



Here 



(A6) 



(A7) 



8ss' 



(i-5 ss ,) 2 + (5 ss ,y 



u, = s 



2 i 2 i a 2 2 , 4 

oo + q + s Acu q + fi 



uj 2 + s'n 2 -(g- uJsf 
Aqu s 
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1/2' 



1/2 



(A8) 



For p 2 < mq, which is the limit taken in the text, Eqs. (A8) reduce to 

uj s ~ sm + q (A9) 

and 

6ss , ~ SS 'V 2 ^ (A10) 

4g(m + sq) I6q 2 (m + sq) 2 

The energy density computed by just substituting Eq. (A6) into the Hamil- 

tonian density is 

p = Po + Re( P 2 q e 2iqt ) + ... (All) 

where the ellipses represent more rapidly varying terms (eg. terms varying 
with frequency u + — U- or 2u s ) and 



Po = 



^0 » 



Piq = p 2 («++ + ia + -) (a_ + + ia — ) . (A12) 

The approximately 'linearly polarized' modes with arg(x) ~ 6/2 correspond 

to (ct_ + + ia )/(a ++ + ia + _) = e* e , whereas the approximately 'circularly 

polarized' modes correspond to the cases (a ++ + ia + ) = and (a_+ + 

ia ) = 0. 

For linearly polarized modes 

p oc (u% + ^/i 2 cos(2gt)^) (A13) 

and the energy per quantum is 

lu 2 

E = uj + -— cos(2gt) (A14) 
2 ujq 

as given in the text. 
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Figure Captions 



Figure 1. Effective coupling to '6- matter' through goldstone boson mixing. 
Figure 2. Mixing of the and Q scalars through a loop. The 77 particle 
carries both a and b charges. 

Figure 3. Similar processes involving the absorption of a) a single axion and 
b) two axions. If the axions come from the coherent state describing the 
classical field oscillation, then the rate for b) exceeds that for a). 
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Fig. 2 
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Fig. 3A 




